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^ ■ Abstract 

In this paper we consider a phase space path integral for general time-dependent quan- 
tum operations, not necessarily unitary. We obtain the path integral for a completely posi- 
tive quantum operation satisfied Lindblad equation (quantum Markovian master equation). 
We consider the path integral for quantum operation with a simple infinitesimal generator. 
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1 Introduction 



Unitary evolution is not the most general type of state change possible for quantum systems. The 
most general state change of a quantum system is a quantum operation [1, 2, 3, 4, 5]. One can 
describe a quantum operation for a quantum system starting from a unitary evolution of some 



closed system if the quantum system is a part of the closed system [6]-[14]. However, situations 

r — ' 

can arise where it is difficult or impossible to find a closed system comprising the given quantum 

j> ! system [15]-[19]. This would render theory of quantum operations a fundamental generalization 

• i—i . 

^ . of the unitary evolution of the closed quantum system. 

The usual models of a quantum computer deal only with unitary quantum operations on pure 
states. In these models it is difficult or impossible to deal formally with measurements, dissi- 
pation, decoherence and noise. It turns out that the restriction to pure states and unitary gates 
is unnecessary [20]. In [20], a model of quantum computations by quantum operations with 
mixed states was constructed. The computations are realized by quantum operations, not neces- 
sarily unitary. Mixed states subjected to general quantum operations could increase efficiency. 
This increase is connected with the increasing number of computational basis elements for the 
Hilbert space. A pure state of n two level quantum systems is an element of the 2 n -dimensional 
functional Hilbert space. A mixed state of the system is an element of the 4 ra -dimensional opera- 
tor Hilbert space. Therefore, the increased efficiency can be formalized in terms of a four-valued 
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logic replacing the conventional two-valued logic. Unitary gates and quantum operations for a 
quantum computer with pure states can be considered as quantum gates of a mixed state quan- 
tum computer. Quantum algorithms on a quantum computer with mixed states are expected to 
run on a smaller network than with pure state implementation. 

The path integral for quantum operations can be useful for the continuous-variable general- 
ization of quantum computations by quantum operations with mixed states. The usual models of 
a quantum computer deal only with discrete variables. Many quantum variables such as position 
and momentum are continuous. The use of continuous-variable quantum computing [21, 22, 23] 
allows information to be encoded and processed much more compactly and efficiently than with 
discrete variable computing. Quantum computation using continuous variables is an alternative 
approach to quantum computations with discrete variables. 

All processes occur in time. It is naturally to consider time dependence for quantum opera- 
tions. In this paper we consider the path integral approach to general time-dependent quantum 
operations. We use the operator space [24]-[38] and superoperators on this space. The path inte- 
gral for unitary evolution from the operator (Liouville) space was derived in [33]. The quantum 
operation is considered as a real completely positive trace-preserving superoperator on the op- 
erator space. We derive a path integral for a completely positive quantum operation satisfied 
Lindblad equation (quantum Markovian master equation) [39, 40, 41, 42, 18, 38]. For example, 
we consider a path integral for a quantum operation with a simple infinitesimal generator. 

In section 2, the requirements for a superoperator to be a generalized quantum operation are 
discussed. In section 3, the general Liouville-von Neumann equation and quantum Markovian 
(Lindblad) master equation are considered. In section 4, we derive path integral for quantum 
operation satisfied Liouville-von Neumann equation. In section 5, we obtain a path integral for 
time-dependent quantum operation with an infinitesimal generator such that the adjoint gener- 
ator is completely dissipative. In section 6, the continuous-variables quantum computation by 
quantum operations with mixed states is discussed. In the appendix, the mathematical back- 
ground (Liouville space, superoperators) is considered. 

2 Quantum operations as superoperators 

Unitary evolution is not the most general type of state change possible for quantum systems. 
The most general state change of a quantum system is a positive trace-preserving map which is 
called a quantum operation. For the concept of quantum operations, see [1, 2, 3, 4, 5]. 

A quantum operation is a superoperator £ which maps the density matrix operator \p) to 
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the density matrix operator \p'). For the concept of superoperators and operator space see the 
appendix and [24]-[38] 

If \p) is a density matrix operator, then £\p) should also be a density matrix operator. Any 
density matrix operator p is a self-adjoint {p\ = p t ), positive (p t > 0) operator with unit trace 
(Trpt = 1). Therefore, the requirements for a superoperator £ to be the quantum operation are 
as follows: 

1. The superoperator £ is a real superoperator, i.e. (<?(A)) = £{A^) for all A. The 
real superoperator £ maps the self-adjoint operator p to the self-adjoint operator £{p): 

(£( P )y = £{p). 

2. The superoperator £ is a positive superoperator, i.e. £ maps positive operators to positive 
operators: £(A 2 ) > for all A ^ or £(p) > 0. 

3. The superoperator £ is a trace-preserving map, i.e. (I\£\p) = (£^(I)\p) = 1 or £\l) = 
I. 

We have to assume the superoperator £ to be not merely positive but completely positive 
[43]. The superoperator £ is a completely positive map of the operator space, if 

n n 
k=l 1=1 

for all operators A k , B k and all n. 

Let the superoperator £ be a convex //near map on the set of density matrix operators, i.e. 

s s 

where all \ s are < \ s < 1 and ^2 S X S = 1. Any convex linear map of density matrix 
operators can be uniquely extended to a //near map on Hermitian operators. Note that any 
linear completely positive superoperator can be represented by 

m m 

£ = Y,L Ak R Al : £{p) = Y.A kP Ai 
k=l k=l 

If this superoperator is a trace-preserving superoperator, then 

m 

E 4A k = L 

k=l 
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The restriction to linear quantum operations is unnecessary. Let us consider a linear real 
completely positive superoperator £ which is not trace-preserving. Let (I\£\p) = Tr(£(p)) be 
the probability that the process represented by the superoperator £ occurs. Since the probability 
is non-negative and never exceed 1, it follows that the superoperator £ is a trace-decreasing 
superoperator: < (I\£\p) < 1 or £\l) < I. In general, any real linear completely positive 
trace-decreasing superoperator is not a quantum operation, since it can be not trace-preserving. 
The quantum operation cannot be defined as a nonlinear trace-preserving operation N by 

ti\p) = £\p){I\£\p)- 1 or ti{p) = -^— (1) 

where £ is a real linear completely positive trace-decreasing superoperator. 

All processes occur in time. It is naturally to consider time dependence for quantum op- 
erations £(t,t ). Let the linear superoperators £(t,t ) form a completely positive quantum 
semigroup [42] such that 

j t £(t,t )=A t £(t,t ), (2) 

where A* is a completely dissipative superoperator [39, 42, 19]. We would like to consider the 
path integral for quantum operations £(t, t ) with infinitesimal generator A, where the adjoint 
superoperator A^ is completely dissipative, i.e. 

A-\A k Ai) - A*(A k )Ai - A k k\A{) > 0. 



3 Evolution equations 

An important property of most open and dissipative quantum systems is the entropy variation. 
Nevertheless the unitary quantum evolution of a mixed state g t described by von Neumann 
equation 

leaves the entropy < S >= —Tr(g t lng t ) unchanged. Therefore, to describe general quantum 
systems, one normally uses [47, 38] a generalization of (3). 

To describe dissipative quantum systems one usually considers [47] the following equation: 

^ = ~[H iet ]+D(e). (4) 
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3.1 Liouville-von Neumann equation 

Let us consider a generalization of equation (3). The Liouville-von Neumann equation [16, 38, 
44] can be represented as the linear equation 

Using the superoperator formalism, this equation can be rewritten in the form 

j t \Q t ) =k t \Q t ). (6) 

The superoperator language allows one to use the analogy with Dirac's notations. This leads 
quite simple to the derivation of the appropriate equations. 

Here A t is a linear Liouville superoperator on the operator space TL. For the Hamiltonian 
(closed) quantum systems (3) this superoperator is defined by the Hamiltonian H: 

A t = ~{L H -R H ). (7) 
For equation (4) the Liouville superoperator has the form 

k t = - l -{L H -R H ) + D. (8) 

In general, the operator \g t ) is an unnormalized density matrix operator, i.e. Trg t = (I\g t ) ^ 1. 
Equation (6) has a formal solution 

\Qt)=£(t,t )\ Qto ), (9) 
where £(t, t ) is a linear quantum operation defined by 

£(Mo) = Texp C drk T . (10) 

J to 

The symbol T is a Dyson's time-ordering operator [45]. The quantum operation (10) satisfies 
the Liouville-von Neumann equation (2). We can define a normalized density matrix operator 

\Pt) by 

i \ i \( T \ \-i i \ £(PM)\q) 

\pt) = \Qt)(I\Qt) , or \ Pt ) = 

(I\S(t,t )\g) 

i.e. p t = g t /Trg t . The evolution equation for the normalized density matrix operator p t can be 
written in the form 

j t \p t ) = A t \p t )-\p t )(I\A t \ Pt ). (11) 

In general, this equation is a nonlinear equation [19]. A formal solution of equation (11) is 
connected with the nonlinear quantum operation (1) by 

\p t )=AT(t,t )\p t0 ). 
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3.2 Quantum Markovian equation 

Lindblad [39] has shown that there exists a one-to-one correspondence between the completely 
positive norm continuous semigroup of superoperators S(t,t ) and superoperator A such that 
the adjoint superoperator A* is completely dissipative. The structural theorem of Lindblad gives 
the most general form of the bounded adjoint completely dissipative Liouville superoperator A. 
The Liouville- von Neumann equation (11) for a completely positive evolution is a quantum 
Markovian master equation (Lindblad equation) [39, 40, 41]: 

^ = ~[H,Pt] + ^ t([V k Pt,Vt] + [V k ,p t Vti). (12) 

fc=i 

This equation in the Liouville space can be written as 

j t \p t ) = Mpt), 

where the Liouville superoperator A is given by 

A = -j{L H -Rh) + ^ E( 2 ^ fc \t - L v t vl - R vl Rv k ). (13) 

The basic assumption is that the general form of a bounded superoperator A, given by the 
Lindblad theorem, is also valid for an unbounded superoperator [42, 46]. Another condition 
imposed on the operators H, V k , v£ is that they are functions of the observables P and Q (with 
[Q, P] = ihl) of the one-dimensional quantum system. Let us consider V k = a k P + b k Q, were 
k — 1,2, and a k , b k are complex numbers, and the Hamiltonian operator H is 



Then with the notation [46]: 



dqq — y Yl ' d PP - O Yl \ b 

Z fc=l,2 Z fc=l,2 

d PQ = -\ Re { a l b k), A = -7m ( J2 a t b k) 

k=l,2 fc=l,2 

equation (12) in the Liouville space can be written as 

j t \pt) = ^L+Lp + ^-L+L Q - (A - p)L P L+\p t ) + (A + P )L Q L + P \p t ) + 
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+d pp L Q L Q \p t ) + d qq LpL P \p t ) - 2d pq LpL Q \p t ), (14) 
where are the multiplication superoperators defined by 

L A = -^{L A - R A ), L A = -(L A + R A ). 

The properties of these superoperators are considered in the appendix. Equation (14) is a super- 
operator form of the well-known phenomenological dissipative model [47, 46]. 

4 Path integral in the general form 

In the coordinate representation the kernel 

Q(q,q',t) = (q,q'\Qt) 
of the density operator \g t ) evolves according to the equation 

Q{q,q',~t) = j dq dq' S(q,q',q ,q' ,t,t )g(q ,q' ,t ). 

The function 

£(q,q',qo,qo,t,t ) = (q,q'\S(t,to)\q ,q' ) (15) 

is a kernel of the linear quantum operation S(t, t ). Let the Liouville superoperator A t be time 
independent, i.e. the quantum operation £(t, t ) is given by 

£(t,t ) = exp{(t-t )A}. (16) 

Proposition 1. Let {£{t, t ),t > t } be a superoperator semigroup on operator space H 

S(t , t ) = /, S(t, t ) = Sit, tjSih, t ), 

where t > t\ > t such that the infinitesimal generator A of this semigroup is defined by (16). 
Then the path integral for kernel ( 15) of the quantum operation £{t, to) has the following form: 

£(?,?', 5b, ?o» Mo) = J VqVq'VpVp' expj^ dt(^[qp - q'p'} + A s (q, q', p, p')) ■ (17) 

This form is the integral over all trajectories in the double phase space with the constraints that 
q(to) = qo, q(t) = q, q'(to) = q' , q'if) = q' and the measure 



Vq = \{dq(t), Vp=H 



dp(t) 
t 2nh 



The symbol As(q, q', p,p') of the Liouville superoperator A is connected with the kernel A(q, q', y, y') 
by 

^s(q,q',P,P f ) = J dydy'A(q,q',y,y') • exp~[(q - y)p - (q' - y')p'}, 

where A(q, q', y, y') = (q, q'\A\y, y') and 

lr % r , 

k{q,q',y,y') = ^ 27rn y n J dpdp'A s (q,q',p,p') ■ exp-[(q - y)p - (q' - y')p'\. 



Proof. 

Let us derive path integral form (17) for quantum operation (16). 

1. Let time interval [t , t] has n + 1 equal parts 

t-t 

T = . 

n + 1 

Using the superoperator semigroup composition rule 

£(t,to) = £(t,t n )£(t n ,t n -i)...£(ti,to), 
where t > t n > t n -x > ... > t x > t , we obtain the following integral representation: 

£(9, 4o, Mo) = / dq n dq' n ...dq 1 dq' 1 £(q,q' \q n ,q' n ,t,t n ) ... £(q 1 ,q' 1 ,q ,q' ,t 1 ,t ). 

This representation can be written in the form 

„ n n+1 

£(<l,<l',<lo,qo,t,to) = J II dq k dq' k \{ £(q k ,q' k ,q k - U q' k _ 1 ,t k ,t k ^ 1 ). 



Here q n+1 = q and q' n+1 = q'. 
2. Let us consider the kernel 

£(q k 

of the quantum operation £(t k , t k -\). If the time interval [t k -i, t k ] is a small, then in the 
coordinate representation we have 

Q(Qk,q' k ,t k ) = J dqk-xdq'^ £(q k , q' k , q k - U <j k _ x , t k , t fe _i)^(? fe _i, 



and 

Q(Qk,q' k ,tk) = (qk,q' k \Qt k ) = (^^(^^-Olft^!) = 



oo 



„=o n! 



= Q(Qk, q' k , t k -i) + r J dqk^dq'^Aiqk, q' k , q k -i, ^_i) ■ g(q k -i, q' k -i, h-i) + 0(r 2 ) = 

= J dqk-idq'^ (% fc -gfc-i)%/Ui-<?fc)+ rA (9fc, q' k , q k -i, q' k -i)+0(T 2 )^) g(q k -i, q' k _ 1: t k -i). 
Therefore, we have 

£(qk, q'k, qk-i, q' k -i,t k , t k -i) = 5(q k -q' k ) +rA(q k , q' k , qk-i, q' k -i) + 0(t 2 ). 



3. Delta-functions can be written in the form 

dp k dp' k i 

J2^ exp h 



S(q k - qk-i)5{q' k -i - 5*) = / J^j^t ex p\\i.qk ~ qk~i)Pk - (q' k ~ q'k-M}- 
Using the relations 



1 i 

(q,q'\p,p') =< q\p><p'W >= j^y^exp-(qp - q'p'), 

1 i 

{p,p'\q,q') =<p\q >< q'\p' >= JK^Z ex P ~ ^(QP-M^ 
we obtain the symbol A. s (qk, q' k ,p k ,p' k ) of the Liouville superoperator by 

Hqk, q'k, Qk-i, q'k-i) = (qk,q'k\Mqk-i,q'k-i) = 

= J dp k dp' k (qk,q'k\MPk,Pk)(Pk,p'k\qk-i,q'k-i) = 

= J dp k dp' k A s (qk,q'k,Pk,p'k) • (qk, q' k \p k ,p' k )(p k ,p' k \q k -i, q'^d = 

1 r i 
= {27ch) 2n J dpkdp ' k As ( qk > q 'kiP^P'k) ■ exp-[{q k - q k -i)Pk ~ {q' k ~ <&-M- 

4. The kernel of the quantum operation S(t k , t k -i) is 

£(qk, q' k , q k -i, q' k -i, t k , tfc-i) = 

= (2nfi) 2n I dpkdp 'k{ l+T ^s{qk,q' k ,p k ,p' k )+0{T 2 )^exp^^ 



5. The kernel of the quantum operation S(t, t ) is 

£(q,Q',qo,Qo,t,t ) = 
f n n+1 dp k dp' i n+1 r 

J II d( lkdq' k II (o^ftjt ex Pj- E - 9fc-i)i>* - («k - 0fc-iK 

n+l 

■n(l-^A 5 (g fc , g ;,p fc ,^)+0(r 2 )). 

k=l 

6. Using 

A/, , ,. A ,A fc , 



lim TTflH -) = lim TT exp( — ), 

n— >oo - L - L 77 ra^oo - L - L n 

k=l U k=l U 

we obtain ^ 

£(Q,Q',Qo,Qo,t,t ) = f IJ n ' ' " ' 



" +1 dp fc dj / A 

/, 1 fcii ( 27r ^) 



2ra 



■exp 



/KQk-qk-i q'k~q'k-i m , a c / 

'2^ r (d Pfe Pk\+ ^s{qk,q k ,Pk,P k ))- 

k=l T 

7. Let q k , q' k , p k , p' k be the values of the functions q(t), q'(t), p(t), p'(t) and t k = t + kr, 
i.e. 

qk = q(t k ), q'k = q'(t k ), p k = p(t k ) , p' k = p(t k ), 

where k — 0, 1, 2, n, n + 1. Using 

lim gfc gfc ~ 1 = gfo-i), lim £ A(f fc )r = / dtA(t), 
we obtain the kernel of the quantum operation S(t, t ) in the path integral form (17). □ 



Corollary. If the dissipative quantum evolution is defined by equation (4), then the path 
integral for the quantum operation kernel has the form 



£(q, q', qo, q' Q , t, t ) = J VqVpVq'Vp' exp(^(A(q,p) - A(p', q')) + D(q, q',p,p') 
Here A(q, p) and A(p', q') are action functionals defined by 

A(q,p) = f dt(qp - H(q,p)), A(p' , q') = f dt(q'p' - H(p', q')). (18) 
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The functional D(q, q', p, p') is a time integral of the symbol Dg of the superoperator D. 

The functional D(q, q', p, p') describes the dissipative part of evolution. 

Corollary. If the quantum system has no dissipation, i.e. the quantum system is a closed 
Hamiltonian system, then D(q, q',p,p') = and the path integral for the quantum operation 
can be separated 

£(q,q',qo,Q'o,t,to) = U*(q,q ,t,t )U(q',q' ,t,t ), 

where 

U(q',q' ,t,t ) = J Vq'Vp' exp-^A(p' , q'), 
U*(q,qo,t,t ) = J VqVpexp^A(p,q). 

The path integral for the dissipative quantum systems and the corresponding quantum oper- 
ations cannot be separated, i.e. this path integral is defined in the double phase space. 

5 Path integral for completely positive quantum operation 

Let us consider the Liouville superoperator (13) for the Lindblad equation. 

Proposition 2. Let {£(t, t ),t > t } be a completely positive semigroup of linear real trace- 
preserving superoperators such that the infinitesimal generators A of this semigroup are defined 
by (13). Then the path integral for the kernel of the completely positive quantum operation 
£(t,t ) has the form 

S(q,q',q ,q' ,t,t ) = J VqVpVq'Vp' F(q,q',p,p') exp^(A(q,p) - A(p' , q')) , (19) 

where A(q,p) and A(p', q') are action functionals (18) and the functional J-{q, q',p,p') is de- 
fined as 

F(q,q',P,p')=exp-— / dt ]T ((V^V k )(q,p) + (V^V k )(p', q') -2V k (q,p)V^p' , q')). (20) 



11 



Proof. The kernel of superoperator (13) is 

A(g, q', y, y') = (q,q k \A\y,y') = ~(< y'\q' >< q\H\y > - <q\y >< y'\H\q' >) + 

-i m -i m 

+tE< >< v'MW > -~ £(< rfW >< Q\vlv k \ y >-< q \ y >< y'Mv k \ q ' >). 

k=l k=l 

The symbol A s (q, q', p, p') of the Liouville superoperator A can be derived by 

A(q,q',y,y') = Jdpdp'(-^-(< y'\p' >< p'\q' >< q^p >< p\y > - 
- < q\p >< p\y >< y'\p' >< p'\H 2 \q' >) + 

+* E < <i\ v k\p><p\y >< y'W ><p'\v k W >), 

k=l 

were the operators Hi and H 2 are defined by the relations 

i m i m 

Hl = H--^V k %, H 2 = H+-^V k %. 

z k=i z k=i 

Then the symbol A s (q, q', p, p') of the Liouville superoperator (13) can be written in the form 

m 

A s (q,q',p,p') = --(H^p) - H 2 (p',q') + i"£V k (q,p)Vi(j/ ,q')), 

k=l 

or 

A s (q,q',p,p') = ~[H(q,p)-H(p',q')}- 

i m 

"2ft E((^*)(?,P) + (V k %)(p',q') - 2V k (q,p)V k \ P ',q')), 

k=l 

where H(q,p) is a gp-symbol of the Hamilton operator H and H(p, q) is a pg-symbol of the 
operator H . 

In the Hamiltonian case (V k = 0), the symbol is given by 

A s (q,q',p,p') = ~[H(q,p) - H(p',q')}. 
The path integral for a completely positive quantum operation kernel has the form 

£(q,q',qo,q'o,t,t ) = J VqVpVqVp F(q,q',p,p') exp^(A(q,p) - A(p',q')). 

Here A(q,p) and A(p',q') are action functionals (18), and the functional jF(g, q',p,p') is de- 
fined by (20). □ 
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The functional T(q, q', p, p') describes the dissipative part of the evolution and can be called 
a (double) phase space influence functional. The completely positive quantum operation is de- 
scribed by the functional (20). 

Corollary. For the phenomenological dissipative model (14) the double phase space path 
integral has the form (19) with the functional 

1 ft / 

J 7 (q,q',p,pf)=exp-z ^ dt[2d qp (q - q')(p - p')- 
-d qq (p ~ p'f - d pp (q - q'f + ih\(pq' - qp') + ihfi(q'p - qpfj . (21) 



Using the well-known connection between the phase space (Hamiltonian) path integral and 
the configuration space (Lagrangian) path integral [51, 52], we can derive the following propo- 
sition. 

Proposition 3. If the symbol A s (q,q',p,p') of the Liouville superoperator can be repre- 
sented in the form 

As(q,q',P,p') = ~[H(q,p) - H(p',q')} + D s (q,q',p,p'), (22) 

where 

H(q,p) = -aki 1 (q)PkPi ~ h(q)p k + c(q), (23) 
D s (q,q',p,p') = -dk(q,q')pk + d' k (q,q')p' k + e(q,q'), (24) 

then the double phase space path integral (17) can be represented as a double configuration 
phase space path integral 

£(q, q', q'o, *, *o) = / VqVq' F(q, q') exp % -(A(q) - A(q')), (25) 

where 

f l 1 1 

A(q)= / dt£(q,q), £(q, q) = -a k i(q)qkqi+aki(q)bk(q)qi+-aki(q)b k (q)bi( y q)-c( y q). (26) 

Jto 2. 2 

This Lagrangian £(q, q) is related to the Hamiltonian (23) by the usual relations 

£{q,q) = qkPk- H(q,p), Pk = ^r- 
13 



Proof. Substituting (22) in (17), we obtain the kernel of the corresponding quantum opera- 
tion. Integrating (17) in p and p', we obtain relation (25) with the functional 

F(q,q') = exp— jf dt(d k (q, q')a k i(q)(bi(q) - —d t (q, q'))+ 
+d' k (q,q')a kl (q')(b l (q') - ^d\(q, q')) + e(q, q') + 5(0) A(q, q')) , 

where 

A(q,q') = - — (ln[det(a kl (q))} + ln[det(a kl (q'))fj . 

□ 

In equation (25) the functional T(q, q') can be considered as the Feynman- Vernon influence 
functional. It is known that this functional can be derived by eliminating the bath degrees of 
freedom, for example by taking a partial trace or by integrating them out. The Feynman- Vernon 
influence functional describes the dissipative dynamics of open systems when we assume the 
von Hove limit for a system-reservoir coupling. One can describe a quantum system starting 
from a unitary evolution of some closed system "system-reservoir" if the quantum system is a 
part of this closed system. However, situations can arise where it is difficult or impossible to 
find a closed system comprising the given quantum system [6]-[14]. 

The Feynman path integral is defined for configuration space. The most general form of 
quantum mechanical path integral is defined for the phase space. The Feynman path inte- 
gral can be derived from the phase space path integral for the special form of the Hamiltonian 
[51, 52, 53, 54, 55]. It is known that the path integral for the configuration space is correct 
[51, 52] only for the Hamiltonian (23). The Feynman- Vernon path integral [6] is defined in the 
double configuration space. Therefore, this path integral is a special form of the double phase 
space path integral (17). The Feynman- Vernon path integral is correct only for the Liouville su- 
peroperator (22), (23), (24). Note that the symbol A s (q, q', p, p') for most of the dissipative and 
non-Hamiltonian systems (with completely positive quantum operations) cannot be represented 
in the form (22). 

Corollary. In the general case, the completely positive quantum operation cannot be rep- 
resented as the double configuration space path integral (25). 
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For example, the double phase space path integral (21) for the phenomenological dissipative 
model (14) has the term pp'. Therefore this model and the Liouville symbol A s (q, q', p, p') for 
this model cannot be represented in the form (22), (23), (24). 

6 On the continuous-variables quantum computation by quan- 
tum operations with mixed states 

The usual models of a quantum computer deal only with the discrete variables, unitary quantum 
operations (gates) and pure states. Many quantum variables such as position and momentum are 
continuous. The use of continuous-variable quantum computing [21, 22, 23] allows information 
to be encoded and processed much more compactly and efficiently than with discrete variable 
computing. Quantum computation using continuous variables is an alternative approach to 
quantum computations with discrete variables. 

In the models with unitary quantum operations on pure states it is difficult or impossible 
to deal formally with measurements, dissipation, decoherence and noise. It turns out that the 
restriction to pure states and unitary gates is unnecessary [20]. In [20], a model of quantum 
computations by quantum operations with mixed states was constructed. It is known that the 
measurement is described by quantum operations. The measurement quantum operations are the 
special case of quantum operations on mixed states. The von Neumann measurement quantum 
operation as a nonlinear quantum gate is realized in [20] . The continuous quantum measurement 
is described by the path integrals [56, 57, 58, 59]. Therefore, the path integral for quantum op- 
erations can be useful for continuous-variables quantum operations on mixed states. Quantum 
computation by quantum operations with mixed states is considered [20] for discrete variables 
only. Some points of the model of the continuous-variable quantum computations with mixed 
states are considered in this section. The double phase space path integral can be useful for the 
continuous-variables quantum gates on mixed states. 

The main steps of the continuous-variables generalization of quantum computations by 
quantum operations with mixed states are following. 

1. The state \p(t)) of the discrete-variable quantum computation with mixed states [20] is a 
superposition of basis elements 

N-l 

\p(t)) = E \l*)Pv(t), (27) 

H=0 



15 



where p^(t) = (p\p(t)) are real numbers (functions). The basis \p) of the discrete- 

— (n) 

variable Liouville space H is defined [20] by 

\fi) = lUi-Hn) = = ^pT^i ® •- ® (28) 

where are Pauli matrices, iV = 4™, each pi £ {0, 1, 2, 3} and 

N-l 
At=0 

is the discrete- variable computational basis. 

The state \p(t)) of the continuous-variable quantum computation at any point of time can 
be considered as a superposition of basis elements 

\p(t)) — dx dx'\x, x')p(x, x', t), (30) 



where p(x,x',t) = (x,x'\p(t)) are the density matrix elements. The basis \x,x') of the 
continuous-variable operator space H is defined by \x, x') = \ \x >< x'\), where 

(x, x'\y, y') — S(x — x')S(y — y'), J dx J dx'\x,x')(x,x'\ — I (31) 

can be considered as a continuous-variable computational basis. 

2. In the discrete-variable computational basis \p) any linear quantum operation £ acting on 
n-qubits mixed (or pure) states can be represented as a quantum four- valued logic gate 
[20]: £ on n-ququats can be given by 

N-l N-l 

£= E E^MH ( 32 ) 

where N = A n , 

= —Tr(a^£(a u j), (33) 

anda M = cr Ml <g> ... <g> a^. 

In the continuous-variable computational basis \x,x') any linear quantum operation £ 
acts on mixed (or pure) states can be represented as a continuous-variable quantum gate: 

£(*2,*i)=y dx dx' dy dy' £(x,x',y,y',t 2 ,t 1 ) \x,x')(y,y'\, (34) 
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where 8(x ) x', y, y\ t 2 , h) = (x, x'\£(t2, h) \y, y') is a kernel of the real trace-preserving 
positive (or completely positive) superoperator £(t 2 , h) This quantum operation can be 
considered as a continuous-variable quantum gate. 

The continuous quantum measurement which is described by the path integral [56, 57, 
58, 59] is the special case of the continuous-variable quantum gate. The path integral for 
the quantum operations can be useful for all continuous-variables quantum operations on 
mixed states. 

3. Many quantum variables, such as position and momentum are continuous. The use of 
continuous-variable quantum computing [21, 22, 23] allows information to be encoded 
and processed much more compactly and efficiently than with discrete variable comput- 
ing. 

Mixed states subjected to the general quantum operations could increase efficiency. This 
increase is connected with the increasing number of computational basis elements for 
operator Hilbert space. A pure state of the quantum systems is an element of functional 
Hilbert space H. A mixed state of the system is an element \p) of the operator Hilbert 
space H. A mixed state of the system can be considered as an element p(x, x', t) of the 
double functional Hilbert space H^H. 

The use of continuous-variable quantum computation by quantum operations with mixed 
states can increase efficiency compared with discrete variable computing. 

7 Conclusion 

The usual quantum computer model deal only with the discrete variables, unitary quantum 
operations and pure states. It is known that many of quantum variables, such as position and 
momentum are continuous. The use of continuous-variable quantum computing [21, 22, 23] 
allows information to be encoded and processed much more efficiently than in discrete-variable 
quantum computer. Quantum computation using continuous variables is an alternative approach 
to quantum computations with discrete variables. 

The quantum computation by quantum operations with mixed states is considered in [20]. 
It is known that the measurement is described by quantum operations. The measurement quan- 
tum operations are the special case of quantum operations on mixed states. The von Neumann 
measurement quantum operation as a nonlinear quantum gate is realized in Ref. [20]. The con- 
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tinuous quantum measurement is described by the path integrals [56, 57, 58, 59]. Therefore, the 
path integral for quantum operations can be useful for continuous-variables quantum operations 
on mixed states. Quantum computation by quantum operations with mixed states is considered 
[20] only for discrete variables. The model of continuous-variable quantum computations with 
mixed states will be suggested in the next publication. The double phase space path integral can 
be useful for continuous- variables quantum gates on mixed states. 

Let us note the second application of double phase space path integral. The path integral 
formulation of the quantum statistical mechanics leads to the powerful simulation scheme [60] 
for the molecular dynamics. In the past few years the statistical mechanics of non-Hamiltonian 
systems was developed for the molecular dynamical simulation purpose [61, 62, 63, 64, 65]. 
The suggested path integral can be useful for the application in the non-Hamiltoniam statistical 
mechanics of quantum [66, 67] and quantum-classical systems [68, 69]. 
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Appendix 

For the concept of Liouville space and superoperators see [24]-[38]. 

A.l. Operator space 

The space of linear operators acting on a Hilbert space H is a complex linear space H. We 
denote an element A of H by a ket-vector \A). The inner product of two elements | A) and \B) 
of H is defined as (A\B) = Tr(A^B). The norm ||A|| = y/(A\A) is the Hilbert-Schmidt norm 
of operator A. A new Hilbert space H with the inner product is called Liouville space attached 
to Ti or the associated Hilbert space, or Hilbert-Schmidt space [24]-[38]. 

The X-representation uses eigenfunctions \x > of the operator X. In general, the operator 
X can be an unbounded operator. This operator can have a continuous spectrum. This leads 
us to consider rigged Hilbert space [48, 49, 19, 38] (Gelfand triplet) B C H = W C B* and 
associated operator space. The rigged operator Hilbert space can be considered as the usual 
rigged Hilbert space for the operator kernels. 
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Let the set {\x >} satisfy the following conditions: 

< x\x' >= 5{x — x') , J dx\x >< x\ = I. 

Then \x, x') = \\x >< x'\) satisfies 

(x : x'\y : y') — 8{x — x')8(y — y') : j dx J dx'\x,x')(x,x'\ = I. 

For an arbitrary element \A) of H, we have 

\A) = J dx J dx'\x,x')(x,x'\A), (35) 

where (x, x'\A) is a kernel of the operator A such that 

(x,x'\A) = Tr((\x >< x'l^A) = Tr(\x' >< x\A) =< x\A\x' >= A(x,x'). 

An operator p of density matrix can be considered as an element \p) of the Liouville (Hilbert- 
Schmidt) space H. Using equation (35), we obtain 

\p) = j dx J dx'\x, x')(x, x'\p) , (36) 

where the trace is represented by 

Trp = (I\p) = J dx (x,x\p) = 1. 

A.2. Superoperators 

Operators, which act on H, are called superoperators and we denote them in general by the hat. 
A superoperator is a map which maps operator to operator. 

For an arbitrary superoperator A on H, which is defined by A\A) = \A(A)), we have 

(x,x'\A\A) = fdy [ dy'(x,x'\A\y,y')(y,y'\A) = [ dy [ dy'A(x, x', y, y')A(y, y'), 



where A(x, x' ', y, y') = (x, x'\A\y, y') is a kernel of the superoperator A. 

Let A be a linear operator in the Hilbert space H. We can define the multiplication super- 
operators L A and R A by the following equations: 

L A \B) = \AB) , R A \B) = \BA). 
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The superoperator kernels can be easy derived. For example, in the basis \x, x') we have 

(x,x'\L A \B) = J dy j dy'(x,x'\L A \y,y')(y,y'\B) = J dy J dy'L A (x,x',y,y , )B(y,y'). 
Using 

(x, x'\AB) =< x\AB\x' >= J dy J dy' < x\A\y >< y\B\y' >< y'\x' >, 
we obtain the kernel of the left multiplication superoperator 

L A (x,x',y,y') =< x\A\y >< x'\y' >= A(x,y)5(x' - y'). 

A superoperator £ t is called the adjoint superoperator for £ if (£\A)\B) = (A\£(B)) for 
all | A) and \B) from TL. For example, if £ = L A R B , then £t = L A tR B u If £ = L A , then 

& = L AU 

Left superoperators L A are defined as Lie and Jordan multiplication by the relations 

L A B = ~z{AB — BA), L A B = -(AB + BA). 

The left superoperator L A and the right superoperator R A are connected by L A = — R A , L A = 
R A . An algebra of the superoperators L A is defined [50] by 

(1) the Lie relations 

L a .b — L A L B — L B L A , 

(2) the Jordan relations 

^(AoB)oC + ^b^c^a + L\L^L% = L\ oB L~c + L BoC L A + L AoC L B , 

^CAoB)oC + ^b^c^a + L A L^L B = L^L AoB + L B L AoC + L A L BoC , 

T + T+ J-T+T+ -L T + T+ - f + f+-i-f+ f + 

^C^AoB ' ^B^AoC ^A^BoC ~ ^AoB^C + ^BoC^A + ^AoC^Bi 

(3) the mixed relations: 

\-B = L A L B — L B L A , L AoB = L\L B + L B L A , 

-^JioB — ^a^b ~^L B L A , L B L A — L^L^ = ^^.b- 
Here we use the notations 

A ■ B = \{AB - BA), AoB = -(AB + BA). 

tit Zi 
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